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Improving palm oil quality through identification 
and mapping of the lipase gene causing oil 
deterioration 
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The oil palm fruit mesocarp contains high lipase activity that increases free fatty acids and 
necessitates post-harvest inactivation by heat treatment of fruit bunches. Even before heat 
treatment the mesocarp lipase activity causes consequential oil losses and requires costly 
measures to limit free fatty acids quantities. Here we demonstrate that elite low-lipase lines 
yield oil with substantially less free fatty acids than standard genotypes, allowing more 
flexibility for post-harvest fruit processing and extended ripening for increased yields. We 
identify the lipase and its gene cosegregates with the low-/high-lipase trait, providing 
breeders a marker to rapidly identify potent elite genitors and introgress the trait into major 
cultivars. Overall, economic gains brought by wide adoption of this material could represent 
up to one billion dollars per year. Expected benefits concern all planters but are likely to 
be highest for African smallholders who would be more able to produce oil that meets 
international quality standards. 
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Oil extracted from the oil palm (Elaeis guineensis) 1 fruit 
mesocarp is the largest source of vegetable oil, with a 
world market share of ~32% (Fig. la,b). The share of 
palm oil production by small, family-owned, estates is 30% 
worldwide and reaches 80% in Nigeria, Africa's largest producer 2 . 
It is expected that the oil palm share of the market will need to 
increase to meet the doubled demand for vegetable oil predicted 
by 2030 (refs 3-6). One major constraint that weighs heavily on 
yield and harvest management is the rapid oil acidification due to 
the high lipase (HL) activity that releases free fatty acids (FFA) in 
the mesocarp of bruised ripe fruits 7-10 . Acidity has a strong 
impact on the quality of commodity oil because an FFA content 
of > 5% is considered to be unfit for human consumption 11 . 

Lipase activity can hydrolyze 30% of the oil content within 
5min 10 and large-scale exploitation of oil palm became feasible 
only after Fickendey 12 found a century ago that the lipase could be 
inactivated by heat treatment of harvested bunches. However, the 
mesocarp lipase can be active prior to heat treatment and release 
FFA. While FFA level of crude palm oil from large estates is 
usually 3-4% 13 , oil from smallholders shows higher values 
especially in Africa where FFA level reaches 9% in Nigeria 14 and 
can exceed 15% 15 . FFAs have to be removed through an energy- 
costly refining process that leads to oil losses 13 , so that most of 
the palm oil produced by African smallholders is not competitive 
on international markets 14 . To limit acidity, the bunches must 
be harvested before they are overripe, heat- sterilized rapidly 
(ideally within 1 day) to inactivate lipase activity and processed 
immediately. However, bunches also need to be harvested when 
they are as ripe as possible to maximize oil yield. These conflicting 
requirements narrow the harvest window down to ~ 1 week and 
oblige planters to harvest every 7-10 days, as oil palm continuously 
produces fruit bunches yearlong (Fig. la). As a result, more than a 
third of the plantation-operating costs are dedicated to harvest 1 . 
Palm oil is currently the least expensive oil on the market, due to 
high productivity and low labour costs in the producing countries, 
mainly Indonesia and Malaysia 4 . This will most likely change in 
the future as oil palm harvest is largely non-mechanized and 
labour-intensive, and wages are steadily increasing as Southeast 
Asia develops 1 ' 4 . It is therefore important that management 
improves to keep harvest-associated costs as low as possible. 

While oil acidification concerns all planters, it is more acute 
for smallholders who may not have access to efficient milling 



facilities close to their plantation and whose overripe bunches 
may be rejected by mills. Indeed, the problem extends from the 
whole industry down to the consumer in those African countries 
where crude, unrefined, palm oil is a major commodity oil. The 
economic importance of palm oil acidity made the mesocarp 
lipase the most widely studied enzyme of oil palm 16 , as early as in 
1910 (ref. 12). However, the source of this lipase activity has not 
yet been identified. 

We previously identified trees with a low mesocarp lipase 
activity . Here we show that the oil extracted from these trees, 
using different procedures that mimic standard processes, have 
significantly lower acidity than that from the standard, HL 
genotypes currently cultivated. We also demonstrate that the 
low- /high -acidity trait is monogenic and identify the lipase 
responsible for oil acidity. Finally, we show that the lipase gene 
cosegregates with the low (LL)/HL trait, making it a useful DNA 
marker for marker- assisted selection of the trait, a key technique 
for fast breeding of crop trees 17 . 



Results 

Oil FFA content is lower in LL lines. To assess the performance 
of LL trees, we used either elite tenera material (Fig. 1) from a 
Deli (Dl) x La Me (LI) cross already used for commercial oil 
production or progenies of its dura Deli (Dl) parent. The oil 
extraction ratio was 22.15% to 21.73% in HL and LL lines from 
Dl, respectively, thus showing no influence of the LL trait on oil 
yield (Table 1). 

When ripe bunches were heat-sterilized a few hours after 
harvest, the oil extracted from standard fruits contained 
0.66 ± 0.36% (s.d.) FFA, while the oil from LL fruits showed less 
than half FFA, down to 0.30 ± 0.02% (Fig. 2a). These ideal 
conditions are rarely met in practice and the average acidity of 
crude palm oil from large estates is usually ~ 3-4% 13 and can 
exceed 10% in small plantations. Therefore, we further focused on 
two critical steps of fruit processing that increase oil acidity: 
delayed harvesting of bunches and post-harvest delayed 
processing. When bunches were left on trees up to 30 days 
after ripeness, the oil from standard trees contained > 12% FFA, 
whereas that from LL trees showed only half that value (Fig. 2b). 
Therefore, it is clear from our data that the use of LL genotypes 
should allow significant extended ripening without a negative 
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Figure 1 | Oil palm is a major oil-crop, (a) Palm tree with several bunches at different stages of ripeness; a palm tree continuously produces bunches 
yearlong (8-28 bunches per year, depending on age). While most cultivated trees are derived from Deli germplasm, several other germplasms such as 
La Me are significantly used by breeders, (b) Open tenera fruits showing fleshy mesocarp from which palm oil is extracted and kernel that contains an 
oil with shorter fatty acids. The amount of mesocarp per fruit is heavily influenced by the Sh gene that governs shell thickness. Tenera trees, which are the 
most productive and most widely cultivated, are heterozygous Sh+/ — . Dura trees are Sh+/+ and contain a thicker shell and a thinner mesocarp, 
thus producing less oil. (c) World vegetable oil production by crop, showing the increasing share of palm oil (FAOSTAT, http://faostat.org, Dec 2011). 
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Table 1 | Yield is not affected in low-lipase lines. 

Line Fruit to bunch Wet mesocarp to fruit Oil to dry mesocarp Mesocarp water content Oil to bunch 

High Lipase 59.27 ±0.58 68.18 ±0.77 79.26 ±0.28 31.62 ±0.31* 21.73 ±0.28 

Low Lipase 59.41 ±0.89 70.12 ±1.18 79.63 ±0.43 33.18 ±0.47* 22.15 ±0.43 

Ratios of main components of yield are in % ± s.e. Means are from 301 bunches from 21 high-lipase trees and 126 bunches from 9 low-lipase trees. A one-way analysis of variance (F test, factor group) 
was carried out after having tested the model s validity (normality, outliers). 

indicates a statistically (P = 0.009) significant difference between low- and high-lipase lines. Only mesocarp water content appears to differ slightly between the two lines. 
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Figure 2 | LL genotype yields oil of lesser acidity than standard HL genotypes, (a) Ripe bunches (1-5 fallen fruits) were heat-sterilized a few hours after 
harvest. Oil was extracted immediately and its acidity titrated, expressed as %FFA). n = 4 for LL fruits and 8 for HL (b) Overripe bunches were left 
up to 30 days on the tree after ripeness (1-5 fallen fruits) before harvest. Bunches and fallen fruits were then sterilized and processed. (c,d) Fruits were 
sprayed with water (control, c) or with fungicide (d), sun-dried and stored up to 12 days before heat sterilization and oil extraction; n = 3-5. Analyses 
of variance were carried out according to the general linear model and differences between standard and LL materials were significant at P< 0.001 
(a-c) and P< 0.012 (d) (Tukey's test). Bunches were harvested from the segregating D1 x L1 cross (a,c and d) or from selfed D1 progenies (LL P4 and 
standard P5, see plant material in Material and Methods section) (b). Error bars represent s.d. 



impact on oil quality. This should increase oil yield as oil content 
of the bunches increases by ~ 7% after first fruit drop starts 18 ' 19 . 

Another critical step is the delay between harvest and bunch 
processing, especially in small estates distant from mills. When 
ripe bunches were processed 3 days after harvest, FFA content of 
oil from the standard genotype was 1.50 ±0.51% compared 
with only 0.35 ± 0.09% from LL trees. The latter value is 
almost identical to that observed when bunches are processed 
immediately after harvest. When processing was further delayed, 
FFA content was still <5% at 12 days in LL lines (Fig. 2c). This 
result clearly indicates that LL lines will allow larger delays for 
processing bunches. This should strongly benefit smallholders 
distant from mills in areas lacking good transportation 
infrastructures 2 , as is frequently the case in Africa and for those 



who traditionally leave the bunches to ferment for several days 
before processing to allow fruit loosening 20 . 

There is a controversy as to the respective contributions of 
endogenous and fungal lipases to oil acidity 21 . Treating bunches 
with fungicide upon harvest (Fig. 2d) led to reduced oil acidity, 
with the oil from LL fruits containing 3% FFA after 12 days. 
Again, the oil from LL trees contained 30% less FFA compared to 
that of standard trees. Our data show that, although both 
mesocarp and fungal lipases are contributors to oil acidity, the 
endogenous lipase is clearly the dominant contributor to the 
immediate post-harvest rise in FFA. To conclude, in conditions 
known to favour high levels of FFA, LL genotypes always showed 
a delayed development of acidity when compared to the standard 
genotype. Therefore, clearly the LL genotypes can produce an oil 
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of significantly lower acidity and therefore higher quality than 
that from presently used genotypes. 

Mesocarp lipase abundance correlates with lipase activity. We 

devised a functional proteomic approach to identify lipases from 
a crude protein extract. The approach is based on the use of an 
inhibitor, tetrahydrolipstatin , that can bind covalently to the 
lipase catalytic serine residue. Thus, a crude protein extract from 
ripe mesocarp from the HL genotype was incubated with 
radiolabeled tetrahydrolipstatin. A single protein band of 
55kDa was subsequently detected by autoradiography (Fig. 3a). 
De novo sequencing (Supplementary Fig. SI A) showed significant 



homology to a protein with high similarity to castor bean acid 
lipase 23 . The MS/MS spectra allowed us to identify from palm 
sequences 24-26 a cDNA named EgLIPl (GenBank accession 
number JX556215), coding for a protein with strong homology to 
castor bean acid lipase (34.6% identity, 51.8% similarity, 
Supplementary Fig. SIB). Antibodies raised against a specific 
peptide detected a 55-kDa protein in HL mesocarp, whereas no 
protein could be detected in the mesocarp from the LL line 
(Fig. 3b). Immunolocalization data (Fig. 3c) showed that the 
putative lipase is localized adjacent to the outer surface of the 
lipid droplets. Taken together, these data suggest that EgLIPl 
encodes the mesocarp lipase. 
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Figure 3 | Isolation and characterization of the oil palm mesocarp lipase EgLIPl. (a) Mesocarp crude protein extract was incubated with radioactive 
tetrahydrolipstatin and loaded on a 10% SDS-PAGE gel (left panel). The gel was dried and autoradiographed (right panel), (b) Western analysis 
of the EgLIPl protein in the fruit mesocarp from LL and HL lines (D1 progenies), (c) Immunolocalization of EgLIPl. Confocal analysis of mesocarp sections 
stained with Nile red (upper panel) to detect lipids in lipid droplets (Ld) 26 labelled with anti-EgLIPI and DAPI (middle panel) or pre-immune serum 
and DAPI as control (lower panel). Yellow arrow indicates anti-EgLIPI signal (yellow spots) at border of Lds in mesocarp cell. Scale bar is 5|im. (d) Lipase 
activity, lipid, EgLIPl protein and transcript profiles in mesocarp during fruit ripening. Fruit and early germinated seedlings were from the PI progeny. 
DAP: days after pollination, gsl, embryo from germinated seed; gs2, embryo from early post-germinated seed. Error bars represent s.d. 
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To obtain more evidence for a role of EgLIPl in the mesocarp, 
we studied its pattern of expression, protein accumulation 
and activity. Lipase activity assays revealed a > 100-fold increase 
during the maturation and ripening stages 26 of fruit development 
(Fig. 3d), which coincided with the onset of lipid accumulation. 
The EgLIPl transcript peak accumulation and the increase 
in the EgLIPl protein (Fig. 3d, Supplementary Fig. S2) occurred 
concomitant to lipase activity. Quantitative RT-PCR analysis 
of transcripts in select palm tissues indicated that the EgLIPl 
gene is almost exclusively expressed in the mesocarp during 
fruit ripening. EgLIPl transcript was not detected in the embryo, 
kernel, flowers, roots, leaves or the vegetative shoot apex 
(Supplementary Fig. S3). In early post-germinated seedlings, the 
EgLIPl transcript was present, but at 34-fold less than 
that observed in the mesocarp (Fig. 3d, Supplementary 
Fig. S3C). To examine the relationships between EgLIPl and 
the LL trait, lipase activity, transcript accumulation and the 
amounts of lipase protein were analyzed in the P2 segregating 
progeny (Fig. 4a, Supplementary Fig. S2, Supplementary Fig. S4). 
Low oil acidity was associated with LL activity, and 
barely detectable amounts of EgLIPl protein and corresponding 
transcript, whereas high levels of acidity correlated with 
high levels of lipase activity, EgLIPl protein and transcripts. 
These results provide strong evidence that oil acidity is 
directly related to the amount of transcript that encodes the 
EgLIPl lipase. 



The acidity trait cosegregates with mesocarp lipase gene. To 

determine the inheritance of the low- acidity trait, we analyzed the 
Dl self-crossed progeny (S). Two groups were observed, one with 
an average FFA content of 11.6% and the other with 48.5% 
(Fig. 4b). The low-acidity group represented 23.8% of the progeny, 
which is consistent with a 3:1 Mendelian segregation and a 
recessive low-acidity trait. We further analyzed three progenies 
derived from two successive Dl self-crosses (Fig. 4b). One con- 
tained only HL individuals (PI), another only LL individuals (P3) 
and the third progeny (P2) showed a 3:1 segregation pattern 
similar to that of S. These results confirmed that the oil acidity 
trait is monogenic with the low acidity being recessive, and that 
Dl is indeed heterozygous for the trait. Notably, similar results 
were also obtained from Dl x LI material (Supplementary Fig. 
S5), which further indicated that the LL trait of Dl and LI trees 
belonged to the same complementation group. 

Analysis of genomic DNA from S and P2 revealed the existence 
of a (TG) 8 (TA) 16 microsatellite marker locus (mEgCIR_LIP03) 
between 744 bp and 792 bp upstream of the EgLIPl initiating 
ATG codon, with two alleles that showed a few base differences. 
A bulk segregant analysis (BSA), carried out in the P2 progeny 
showed a cosegregation of mEgCIR_LIP03 with the oil acidity 
trait. Two marker alleles (267 bp and 269 bp) were present in the 
bulk for HL activity and one single marker allele (267 bp) in the 
bulk for LL activity (Fig. 4c). Our BSA with 10 genotypes per bulk 
showed a very low probability (2 x 10 ~ 6 ) that this marker locus 
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Figure 4 | The low-acidity trait is monogenic and cosegregates with EgLIPl lipase gene, (a) Lipase activity and EgLIPl protein in P2 progeny. Error bars 
represent s.d. (b) Acidity measurements from the S, PI, P2 and P3 progenies. Selfing Dl led to a progeny of 21 trees (S), 3 of which were further self-crossed 
to give rise to progenies PI, P2 and P3 (see plant material in Material and Methods section). For each progeny, width of Boxplots is proportional to the 
number of trees per class. S and P2 showed a 3:1 ratio consistent with a monogenic recessive trait (S: = 0.016, P = 0.90, d.f. = 1; P2: % 2 = 0.205, 
P = 0.65, d.f. = 1) (c) The low-acidity trait is linked to the EgLIPl locus. The microsatellite marker EgCIR_LIP03 is located in the promoter region of EgLIPl 
and shows two alleles that are visualized on electrophoregrams as 267 and 269-bp-long amplicons. Bulk segregation analysis was carried out on P2 
progeny. Each genomic DNA pool (HL or LL) was represented by 10 randomly chosen trees in each acidity class. 



NATURE COMMUNICATIONS | 4:2160 | DOI: 10.1038/ncomms3160 | www.nature.com/naturecommunications 

© 2013 Macmillan Publishers Limited. All rights reserved. 



5 



ARTICLE 



NATURE COMMUNICATIONS | DPI: 10.1038/ncomms3160 



is polymorphic without a physical link with the major gene 
controlling the oil acidity. The cosegregation provides further 
evidence that EgLIPl is the gene that controls oil acidity in the 
pulp of mature fruits, or is at least closely linked to another gene 
controlling the activity of EgLIPl itself. In either case, the 
mEgCirLIP03 marker or EgLIPl itself provide useful markers for 
introgression of the LL trait into other elite oil palm genotypes. 

Discussion 

Our results indicate the EgLIPl transcript is almost exclusively 
present in the mesocarp, suggesting that EgLIPl function is 
restricted to this tissue (Supplementary Fig. S3). The oil contained 
in the mesocarp of oil palm fruit serves to attract disseminating 
animals and does not participate in the metabolism of other 
tissues and organs of the plant. One can speculate that the lipase 
function may enhance the palatability of fat, given this is due in 
part to receptors to fatty acids and not to triacylglycerol 27 , 
therefore increasing fruit attractiveness towards disseminating 
animals. As the main components of yield are not affected in LL 
trees (Table 1), the LL allele was evidently not selected against by 
breeders. As Deli genotypes originated from a narrow (4 trees) 
genetic basis, the LL allele is likely present in a significant number 
of elite trees currently used for seed production. Indeed, a 
previous study allowed us to identify 4 LL trees out of 37 (ref. 11). 
We are currently screening a wide range of material for the LL 
allele to assess whether the microsatellite marker can be used in 
genotypes other than Deli and to increase our mapping 
population to achieve fine mapping of the trait. Preliminary 
mapping data confirm BSA results (N. Billotte, unpublished data). 

The availability of a molecular marker for the trait will allow 
breeders to rapidly identify potent genitors that carry a LL allele, 
and to introgress the trait into main elite genotypes during the 
next breeding cycle for seed production. The low-lipase genotypes 
should allow more flexibility in the management of harvest and 
bunch processing. Delaying harvest by 10 days in large estates 
should still reduce losses in the form of FFA, provide an 
additional 2-3% oil yield due to extended ripening and reduce 
harvesting frequency by one-third, thus lowering labour costs. 
Taken together, these impacts should represent an economic gain 
of ~ 1 billion dollars per year, based on today's production. 
Substantially higher gains are anticipated provided the LL trait is 
associated to non- or low-shedding fruits 28 . However, the most 
important benefits are likely to be for African smallholders whose 
improved oil will gain access to international markets. 

Methods 

Plant material. All plant material were derived either from a cross between a tree 
(Dl) from the Deli population and a tree (LI) from the La Me population or from 
progenies of Dl (Supplementary Fig. S6). The offspring of Dl x LI cross 
corresponds to elite material planted worldwide. This cross was used to assess the 
performances of LL trees, unless otherwise indicated. Analyses were conducted 
using progenies from Dl. S progeny resulted from self-crossing of Dl and contains 
21 trees, 5 of which were further self-crossed to give rise to progenies PI (15 trees), 
P2 (104 trees), P3 (17 trees), P4 (75 trees) and P5 (75 trees). These trees are 
cultivated both at CRA-PP (Pobe, Benin) and CEREPAH (La Dibamba, 
Cameroon). A bunch was considered to be ripe when one to five fruits detached. 
Each experiment was performed 3-5 times. 



Identification of LL and HL trees. The amount of endogenous lipase was 
estimated by measuring the percentage of FFA released from a mesocarp 
homogenate for 1 hour after tissue disruption 10 . Only attached, undamaged fruits 
that showed no visual signs of fungal infestation were used. Mesocarp (10 g) was 
ground with mortar and pestle and stored at room temperature for 1 hour. 
The mixture was extracted with hexane (lOmlg -1 of pulp) and debris filtered 
away. Hexane was evaporated and the oil dried in a desiccator for 1 hour. Oil 
acidity was titrated with 0.1 N KOH according to standard procedures 11 . FFA 
content was calculated based on the molecular mass of palmitic acid. 



Oil extraction from bunches. Oil was extracted using procedures similar to those 
used at mills. Bunches were subdivided into spikelets, which were boiled for 1-2 h; 
the fruits detached manually and the oil was extracted using a screw-press. After 
decantation, the oil was skimmed, dried, and its acidity was determined by titration. 
Data were Box-Cox transformed to achieve a normal distribution, ascertained by 
Anderson-Darling test. Analyses of variance were performed according to the 
general linear model. Multiple mean comparisons were made by Tukey's test. 



Fungicide treatment. After harvest, each bunch was divided; one half was sprayed 
with a fungicide solution of DITHANE M45 (mancozeb 80%) at 2gl~ 1 and the 
other was used as a control, and sprayed with water. Both samples were sun-dried 
and kept in a closed room at room temperature for 6, 9 and 12 days before the oil 
was extracted. Multiple mean comparisons were made by Tukey's test. 



Mesocarp constituent analyses. For each stage of development, three independent 
bunches were collected from three individuals of the same progeny. Fifteen 
spikelets were collected in the centre of each bunch, from which five spikelets were 
randomly sampled. Three undamaged fruits were then randomly sampled from each 
of the five spikelets. Their mesocarp was dissected and flash-frozen in liquid 
nitrogen. Total lipids were extracted from 1-g samples of freeze-dried powder using a 
modified Folch method as described previously 29 . 

For all electrophoresis gel analyses, total mesocarp proteins were extracted using 
aqueous phenol 30 . Labelling the lipase with radioactive tetrahydrolip statin was 
carried out as described 30 , except that THL was dissolved in hexane and added to 
the hexane extract of mesocarp without using detergent. The gel was dried and 
submitted to autoradiography. The labelled band was submitted to trypsin 
digestion and analyzed by LC-MS/MS at Bordeaux Proteomic Facility as described 
previously 30 ' 31 . Lipase activity was assayed by continuous titration of released fatty 
acids with a pHSTAT 8 . Antibodies were raised by Biotem France against the 
peptides RDYVNGGRLAKITGKYS and DRKKGFPKELSRND. The peptides were 
conjugated to a protein carrier and injected with Freund complete adjuvant. 
Rabbits were boosted twice and exsanguinated at day 45. 

Total RNA from oil palm tissues was extracted as previously described 32 ' 33 and 
1 ug was used to synthesize cDNA using a first- strand cDNA synthesis kit 
(ImProm-II Reverse Transcription System, Promega). qPCR was conducted on a 
LightCycler 480 (Roche) in 96-well plates in a volume of 10 ul containing 2 ul of 
cDNA diluted 1/100, 1.5 ul of each primer (2uM; 5' - ACTACGTCAACGGTG 
GAAGG-3', S'-CCGTTCCAAATCCAGTTCTT-S') and 5ul SYBR Green 
Mastermix (Roche). PCR was initiated by denaturation at 95 °C for lOmin, 
followed by 45 cycles of 95 °C for 15 s, 60 °C for 15 s and a final extension at 70 °C 
for 1 min. Transcript levels were normalized to that of EgEfal (GenBank accession 
number: AY550990, oligonucleotides 5'-TATCAAAGGATGGGCAGACC-3' and 
5' -TCATCATACCTTGCCTTGG A- 3' ) mRNA, and relative mRNA abundance 
was determined as previously reported 34 . Control using RNA matrices were also 
conducted to validate the absence of DNA in each sample. Each time point was 
replicated three times from 3-5 independent biological samples, and all amplified 
cDNA fragments were sequenced by Beckman-Cogenics. Gene abundance is 
expressed as mean and standard error bars are calculated from three technical 
replicates. 



Immunolocalization of EgLIPl lipase. Mesocarp tissue samples were fixed in 
0.1 M PBS buffer containing 8% (w/v) paraformaldehyde and 4% (w/v) 
^-(S-Dimethylaminopropy^-Ar'-ethylcarbodiimide hydrochloride for 16 h at 4°C. 
Samples were subsequently washed in 0.1 M PBS buffer and cut into 150 um 
sections using a vibratome. Sections were saturated in 0.1 M PBS containing 3% 
(w/v) foetal bovine serum for 24 h at 4°C and transferred into 0.1 M PBS buffer 
containing either anti-EgLIPl antibodies or pre-immune serum (diluted 1/100) for 
24 h at 4 °C. Sections were then washed 4 times in 0.1 M PBS and transferred into a 
0.1 M PBS solution containing goat anti-rabbit IgG antibody Alexa Fluor 545 
conjugated (dilution 1/500) for 3 h. Sections were washed four times in PBS (0.1 M) 
and transferred to PBS buffer containing DAPI (14.3 uM) to visualize nucleus and 
membrane of oil droplets. Observation was performed using a confocal microscope 
(Zeiss LM510 meta) at magnifications x 40/1.2 numeric aperture, at the Cirad 
imagery facility (Montpellier, France). Excitation was performed at 543 nm for 
Alexa Fluor 545 and at 405 nm for DAPI. 



Cloning and sequencing of the EgLIPl genomic region. A bacterial artificial 
chromosome (BAC) library (92,160 clones, ~6 genome equivalents) made from 
the oil palm genotype LI was screened using EgLIPl cDNA as a probe. Three 
positive BAC clones were sequenced (with at least 50-fold coverage) using the 
Roche/454 GS FLX sequencing technique and assembled by Beckman-Cogenics 
(USA). The BAC contigs were assembled into one single haplotype sequence using 
CodonCodeAligner version 3.7.1 (CodonCode Corporation, USA), although the 
tree used as source of DNA proved to be heterozygous for the lipase trait. A search 
of microsatellite loci in the EgLIPl genomic region and the design of PCR primers 
for microsatellite locus mEgCirLIP03 were performed using the SSR Analysis Tool 
(SAT) software 35 . 
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Bulk segregant analysis. Genomic DNA was extracted using the commercial kit 
NucleoSpin 96 Plant according to the manufacturer instructions (Macherey-Nagel, 
Germany). Legitimacy of P2 genotypes was confirmed by ADN ld (Montpellier, 
France) company using 14 E. guineensis microsatellite markers 36 . A BSA was 
performed on the progeny P2 showing a 3:1 Mendelian segregation pattern for the 
oil acidity trait. The mEgCIR_LIP03 marker locus was amplified in the parent of P2 
and in two equimolar mixtures of 10 DNAs each and made from P2 individuals 
randomly chosen within the two extreme phenotypic classes of oil acidity: (i) of HL 
phenotypes (oil acidity > 40%) supposed to be homozygous ( + + ) or 

heterozygous ( H ) for EgLIPl, and (ii) of LL genotypes (acidity <20%) 

supposed homozygous ( ) for that gene. The DNAs were PCR amplified in 

condition of fluorescence 36 ' 37 using the mEgCirLIP03 primers (forward: 5'-TCAAA 
GAAACTTGTAGCATATATCAAA-3'; reverse: 5'-CATCCAGTAAGCTAACA 
CACAAATG-3') in a Mastercycler ep 384 thermocycler (Eppendorf, Germany) on 
25 ng of DNA in a 10 ul final volume of buffer (10 mM Tris-HCl (pH 8), 200 uM of 
dNTP, 0.5 uM MgCl 2 , 0.05% wlv gelatin) containing 0.08 uM of the Ml 3 -labelled 
primer, 0.1 uM of the other primer and 1 U of Taq DNA polymerase. The PCR 
programme was: initial denaturation at 94 °C for 4min; 10 «Touchdown» cycles 
each of 94 °C for 45 s, 55 °C ( - 0.5 °C per cycle) for 1 min, 72 °C for 1 min; and 
then 25 cycles of 94 °C for 45 s, 50 °C for 60 s and 72 °C for 60 s; and a final 
elongation step at 72 °C for 30 min. The PCR products were revealed by capillary 
electrophoresis on an Applied Biosystems 3500xL Genetic Analyzer (Life 
Technologies, USA) and the chromatogram allelic patterns were analyzed under 
the GeneMapper Software version 4.1 (Life Technologies, USA). 



Determination of yield components. Ripe bunches were collected and analyzed 
as previously described 38 . Oil was extracted and quantified as described above. 
Results were analyzed using SAS v9.2 software (MIXED Procedure, SAS Institute, 
http://www.sas.com): a one-way analysis of variance (factor group) was carried out 
and the model's residuals were studied to test the model's validity (normality, 
outliers). The tree/bunch effect was tested with a likelihood ratio test. If significant, 
it was kept in the model as a random effect. The results were expressed using a 95% 
confidence interval and a margin of error of 5%. 

References 

1. Corley, R. H. V. & Tinker, P. B. The Oil Palm 4th edn (Blackwell Science, 2003). 

2. Vermeulen, S. & Goad, N. In Natural Resource Issues Series No. 5 (International 
Institute for Environment and Development, 2006). 

3. Carlsson, A. S., Yilmaz, J. L., Green, A. G., Stymne, S. & Hofvander, P. 
Replacing fossil oil with fresh oil— with what and what for? Eur. J. Lipid Sci. 
Tech. 113, 812-831 (2011). 

4. Carter, G, Finley, W., Fry, J., Jackson, D. & Willis, L. Palm oil markets and 
future supply. Eur. }. Lipid Sci. Tech. 109, 307-314 (2007). 

5. Corley, R. H. V. How much palm oil do we need? Environ. Sci. Policy 12, 
134-139 (2009). 

6. Murphy, D. J. Oil palm: future prospects for yield and quality improvements. 
Lipid Technology 21, 257-260 (2009). 

7. Henderson, J. & Osborne, D. J. Lipase activity in ripening and mature fruit of 
the oil palm— stability in vivo and in vitro. Phytochemistry 30, 1073-1078 
(1991). 

8. Ngando-Ebongue, G. F., Dhouib, R., Carriere, F., Amvam Zollo, P. H. & 
Arondel, V. Assaying lipase activity from oil palm fruit {Elaeis guineensis Jacq.) 
mesocarp. Plant Physiol. Biochem. 44, 611-617 (2006). 

9. Sambanthamurthi, R., Oo, K. C. & Parman, S. H. Factors affecting lipase 
activity in Elaeis guineensis mesocarp. Plant Physiol. Biochem. 33, 353-359 
(1995). 

10. Desassis, A. Palm oil acidification. Oleagineux 12, 525-534 (1957). 

11. Ngando Ebongue, G. F. et al. Identification of oil palm breeding genotypes 
producing an oil with low acid value. Eur. J. Lipid Sci. Tech. 110, 506-509 
(2008). 

12. Fickendey, E. Uber die Verwendungsfahigkeit des Palmols als Speisefett. 
Tropenpflanzer 14, 566-570 (1910). 

13. Gibon, V., De Greyt, W. & Kellens, M. Palm oil refining. Eur. J. Lipid Sci. Tech. 
109, 315-335 (2007). 

14. Ohimain, E. I., Daokoru-Okukole, C, Isah, S. C. & Alaka, E. E. Assement of the 
quality of crude palm oil produced by smallholder processors in rivers state, 
Nigeria. Niger. J. Agri. Food Environ. 8, 28-34 (2012). 

15. Osei-Amponsah, C. et al. Processing practices of small-scale palm oil producers 
in the Kwaebibirem district, Ghana: a diagnostic study. NJAS- Wagen J. Life. Sci. 
60-63, 49-56 (2012). 

16. Sambanthamurthi, R., Sundram, K. & Tan, Y. A. Chemistry and biochemistry 
of palm oil. Prog. Lipid Res. 39, 507-558 (2000). 

17. Mayes, S. et al. In Genomics of Tropical Crop Plants, (eds Moore, P. H. & 
Ming, R.) 371-404 (Springer, 2008). 

18. Rao, V., Gomez, M., Chayawat, N., Chia, C. C. & Lim, C. C. In Proceedings of 
the PIPOC 2001 International Palm Oil Congress (Agriculture) 167-191. 

19. Corley, R. H. V. & Law, I. H. Ripening, Harvesting and Oil Extraction. The 
Planter 77, 507-524 (2001). 



20. Zu, K. S. A., Adjei-Nsiah, S. & Bani, R. J. Effect of processing equipment and 
duration of storage of palm fruits on palm oil yield and quality in the 
Kwaebibrem District, Ghana. Agric. Res. Rev. 1, 18-25 (2012). 

21. Tombs, M. P. & Stubbs, J. M. The absence of endogenous lipase from oil palm 
mesocarp. /. Sci. Food. Agr. 33, 892-897 (1982). 

22. Hadvary, P., Lengsfeld, H. & Wolfer, H. Inhibition of pancreatic lipase 

in vitro by the covalent inhibitor tetrahydrolipstatin. Biochem. J. 256, 357-361 
(1988). 

23. Eastmond, P. J. Cloning and characterization of the acid lipase from castor 
beans. /. Biol. Chem. 279, 45540-45545 (2004). 

24. Bourgis, F. et al. Comparative transcriptome and metabolite analysis of oil palm 
and date palm mesocarp that differ dramatically in carbon partitioning. Proc. 
Natl Acad. Sci. USA 108, 12527-12532 (2011). 

25. Nurniwalis, A. W., Suhaimi, N., Siti Nor Akmar, A., Aminah, S. & Mohamad 
Arif, M. A. Gene discovery via expressed sequence tags from oil palm (Elaeis 
guineensis Jacq.). /. Oil Palm Res. 20, 87-96 (2008). 

26. Tranbarger, T. J. et al. Regulatory mechanisms underlying oil palm fruit 
mesocarp maturation, ripening, and functional specialization in lipid and 
carotenoid metabolism. Plant Physiol. 156, 564-584 (2011). 

27. Khan, N. A. & Besnard, P. Oro-sensory perception of dietary lipids: new 
insights into the fat taste transduction. Biochim. Biophys. Acta 1791, 149-155 
(2009). 

28. Corley, R. H. V. Progress with non-shedding oil palms. Planter 82, 799-807 
(2006). 

29. Laffargue, A., de Kochko, A. & Dussert, S. Development of solid-phase 
extraction and methylation procedures to analyse free fatty acids in lipid- rich 
seeds. Plant Physiol. Biochem. 45, 250-257 (2007). 

30. Dhouib, R. et al. Identification of a putative triacylglycerol lipase from papaya 
latex by functional proteomics. FEBS f. 278, 97-110 (2011). 

31. Deytieux, C. et al. Proteome analysis of grape skins during ripening. /. Exp. Bot. 
58, 1851-1862 (2007). 

32. Argout, X. et al. Towards the understanding of the cocoa transcriptome: 
production and analysis of an exhaustive dataset of ESTs of Theobroma cacao L. 
generated from various tissues and under various conditions. BMC Genomics 9, 
512 (2008). 

33. Morcillo, F. et al. Somaclonal variation in micropropagated oil palm. 
Characterization of two novel genes with enhanced expression in epigenetically 
abnormal cell lines and in response to auxin. Tree Physiol. 26, 585-594 (2006). 

34. Pfaffl, M. W. A new mathematical model for relative quantification in real-time 
RT-PCR. Nucleic Acids Res. 29, e45 (2001). 

35. Dereeper, A., Argout, X., Billot, C, Rami, J. F. & Ruiz, M. SAT, a flexible and 
optimized Web application for SSR marker development. BMC Bioinformatics 
8, 465 (2007). 

36. Roy, R., Steffens, D. L., Gartside, B., Jang, G. Y. & Brumbaugh, J. A. Producing 
STR locus patterns from bloodstains and other forensic samples using an 
infrared fluorescent automated DNA sequencer. /. Forensic Sci. 41, 418-424 
(1996). 

37. Billotte, N. et al. Microsatellite-based high density linkage map in oil palm 
(Elaeis guineensis Jacq.). Theor. Appl. Genet. 110, 754-765 (2005). 

38. Blaak, G., Spamaaij, L. D. & Mendez, T. Breeding and inheritance in oil palm 
(E. guineensis Jacq.). II. Methods of bunch quality analysis. /. W. Afr. Inst. Oil 
Palm Res. 4, 145-146 (1963). 



Acknowledgements 

This study was supported in part by PalmElit SAS, the Aquitaine region and by recurrent 
funding from IRD, Cirad, INRAB, IRAD, CNRS and University of Bordeaux. H.D. and 
G.-F.N.-E. were recipients of SCAC French embassy studentships, G.-F.N.-E. of an 
EMBO fellowship and F.B. of a post-doctoral fellowship from the Aquitaine region. We 
thank G. Piombo for helpful technical assistance and advice on oil acidity measurements 
and V. Poulain for statistical analysis of data. We especially thank S. Dussert for the lipid 
content determination, and A. Crabos and G. Bourdel for their help in RNA analysis. 
Radioactive tetrahydrolipstatin was obtained from Drs P. Hadvary and H. Lengsfeld 
(Roche) and the perl script for de novo sequencing from B. Valot (INRA, Ferme du 
Moulon). We are greatly indebted to J. Ohlrogge for helpful comments on the 
manuscript. 

Author contributions 

F.M., D.C., N.B., G.-F.N.-E., B.N. and V.A. designed the research. F.M., T.J.T., D.C., H.D., 
F.B., R.D., G.-F.N.-E., S.C., T.C., S.E., N.B., M.P. and V.A. performed research. V.A., 
F.M., D.C., N.B. and T.J.T. analyzed the data and wrote the paper. All authors discussed 
the results, and read and approved the final submitted manuscript. 

Additional information 

Supplementary Information accompanies this paper at http://www.nature.com/ 
naturecommunications 



NATURE COMMUNICATIONS | 4:2160 | DOI: 10.1038/ncomms3160 | www.nature.com/naturecommunications 

© 2013 Macmillan Publishers Limited. All rights reserved. 



7 



ARTICLE 



NATURE COMMUNICATIONS | DOI: 10.1038/ncomms3160 



Competing financial interests: The study was funded in part by PalmElit, a 
breeding company that might market low-lipase lines of oil palm in the future. 
Bruno Nouy is a PalmElit employee. All other authors declare no competing financial 
interests. 

Reprints and permission information is available online at http://npg.nature.com/ 
reprintsandpermissions/ 



How to cite this article: Morcillo, F. et al. Improving palm oil quality through 
identification and mapping of the lipase gene causing oil deterioration. Nat. Commu 
4:2160 doi: 10.1038/ncomms3160 (2013). 



This work is licensed under a Creative Commons Attribution- 
NonCommercial-ShareAlike 3.0 Unported License. To view a copy of 
this license, visit http://creativecommons.Org/licenses/by-nc-sa/3.0/ 



8 



NATURE COMMUNICATIONS | 4:2160 | DOI: 10.1038/ncomms3160 | www.nature.com/naturecommunications 
© 2013 Macmillan Publishers Limited. All rights reserved. 



